arXiv:1508.06057v2 [hep-ph] 11 Jan 2016 


Inhomogeneous condensation in effective models for QCD using the finite-mode 

approach 

Achim Heinz“, Francesco Giacosa“’^, Marc Wagner®, Dirk H. Rischke® 

“ Goethe-Universitdt Frankfurt am Main, Institut fur Theoretische Physik, 
Max-von-Laue-Strafle 1, D-60438 Frankfurt am Main, Germany and 
^ Institute of Physics, Jan Kochanowski University, ul. Swietokrzyska 15, 25-406 Kielce, Poland 

We use a numerical method, the finite-mode approach, to study inhomogeneous condensation in 
effective models for QCD in a general framework. Former limitations of considering a specific ansatz 
for the spatial dependence of the condensate are overcome. Different error sources are analyzed and 
strategies to minimize or eliminate them are outlined. The analytically known results for 1 -I- 1 
dimensional models (such as the Gross-Neveu model and extensions of it) are correctly reproduced 
using the finite-mode approach. Moreover, the NJL model in 3 -I- 1 dimensions is investigated and 
its phase diagram is determined with particular focus on the inhomogeneous phase at high density. 

PACS numbers: 12.39.Fe, ll.30.Rd, 11.30.Qc, 12.38.Lg. 

Keywords: nonzero temperature, nonzero density, chiral restoration, inhomogeneous condensation, Gross- 
Neveu model, NJL model. 


I. INTRODUCTION 

Quantum Chromodynamics (QCD) cannot be solved analytically at low energies. However, several aspects of 
QCD can be understood by using effective models which exhibit the same symmetries as QCD, most notably chiral 
symmetry. Some models utilize exclusively hadronic degrees of freedom [such as chiral cr models [H-Q], while others 
feature constituent quarks [such as the Nambu-Jona-Lasinio (NJL) model [MU and the Gross-Neveu (GN) model 
p^]. A model with both hadronic and quark degrees of freedom has also been discussed [3 [ 3 . 

All these effective descriptions of QCD include the spontaneous breaking of chiral symmetry, which implies the 
emergence of a chiral condensate at low temperatures and densities denoted as a. This quantity is represented by 
a nonzero expectation value of a scalar-isoscalar mesonic field in hadronic chiral models or, equivalently, by the 
quark-antiquark expectation value (^/’'0) in quark-based models. 

The chiral condensate is, in general, a function of space, cr = a{x). In principle, the determination of (j{x) is 
straightforward: one has to find the field configuration which minimizes the effective action at a given temperature 
and density. In practice, this task is, however, very difficult. This is why cr is often assumed to be spatially constant. 
This assumption is usually valid in the vacuum and at low densities, but not anymore at high densities. One of 
the simplest non-constant field configurations is the so-called chiral density wave (CDW) which, in chiral hadronic 
models, corresponds to a one-dimensional condensate of the form cr(a; 3 ) = (j)cos{pxz) together with pion condensation, 
TT^ixs) = (j)sin{px 3 ). Various studies have found that the CDW is favorable compared to a constant condensate at 
sufficiently hi gh d ensities [l5l - [3| . Interestingly, a CDW has recently also been obtained within the extended Linear 
Sigma Model [13, which is a general chiral hadronic model with (axial-)vector degrees of freedom 0,i. Moreover, 
inhomogeneous phases were also investigated in Refs. [28l - l3| in the framework of the NJL model as well as in the 
quark-meson model and the skyrmion model. 

A general method to determine space-dependent condensates at non-zero temperature and density has not yet been 
established. There are a few exploratory studies of such methods in the context of the 1-1-1 dimensional GN model, 
using either a lattice regularization |42j [ of the effective action or an expansion in terms of plane waves or hat-like 
localized basis functions [H, 3 , but quite often one uses a specific Ansatz [3l|, 3 ■ More evolved models including 
three spatial dimensions and two-dimensional variations of the condensates and corresponding general methods, 
which are also based on expansions of fermionic fields and condensates in terms of plane waves, have been discussed 
theoretically in Refs. 35 0,3' lo practice, however, due to limited computational resources, investigations have 
again been limited to specific Ansatze, where only a selected set of Fourier modes is considered, e.g. variations in only 
a single spatial dimension. In this respect models for which analytic inhomogeneous solutions are known are extremely 
interesting. This is the case for the 1-1-1 dimensional GN model where a soliton-like solution for the spatial 

dependence of the condensate is found, which is mathematically represented by a Jacobi elliptic function 3 jE 3I' 
Further interesting 1-1-1 dimensional models for which inhomogeneous phases have analytically been determined are 
extensions of the GN model: the chiral Gross-Neveu (yCN) model 35 0|i which has a continuous chiral symmetry, 
and the two-flavor NJL 2 model 3l' These 1-1-1 dimensional models are relevant, because at high densities QCD 
effectively reduces from 3 -I- 1 to 1 -I- 1 dimensions 3) 101H3 ■ 

Thus, while the existence of inhomogeneous phases has been verified by several different approaches, it is highly 
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desirable to develop a general and reliable numerical method to study inhomogeneous condensation, which does not 
require a specific ansatz for the spatial dependence of the condensate. This is the aim of the present work. We 
adapt and extend techniques introduced and explored in Refs. [13,13. We first test the validity and reliability of the 
resulting method, the finite-mode approach, by applying it to 1 -b 1 dimensional models, the GN, the yGN, and the 
NJL 2 model. We correctly reproduce both soliton-like and CDW modulations without supplying any specific ansatz. 

Then we apply the finite-mode approach to study the phase structure of the 3-1-1 dimensional NJL model. Recent 
findings concerning one-dimensional modulations are confirmed. In addition, we determine the shape of the 
so-called inhomogeneous “continent” at high density [H, [H^: in agreement with these works, the phase boundary 
between chirally restored and inhomogeneous phase first increases with temperature. However, for larger chemical 
potential /i it decreases. Thus, the inhomogeneous phase exhibits a shape which is surprisingly similar to that of the 
crystal phase of the GN model. 

The paper is organized as follows. In Sec. |ll] quark-based effective models for QGD in 1 -b 1 and 3 -b 1 dimensions are 
introduced. In Sec. mi Sec. mi Sec. El and Sec. ED the phase diagrams of these models are investigated numerically 
using the finite-mode approach, with particular focus on inhomogeneous condensation. Finally, we present conclusions 
and an outlook in Sec. Em 


II. QUARK-BASED EFFECTIVE MODELS 

In this section we introduce the Lagrangians of the models that we use to investigate inhomogeneous condensation. 
We start with 1 -b 1 dimensional models and then turn to the 3 -b 1 dimensional NJL model. 


A. 1 -b 1 dimensions: the GN model and its extensions 

GN model 

The GN model [H, IH, is a fermionic model that contains only a single quark flavor. In the large-TV 

limit (where N is the number of colors) it exhibits QCD-like features such as asymptotic freedom, dynamical chiral 
symmetry breaking and its restoration, dimensional transmutation, and meson and baryon bound states [bD - l^ . The 
Lagrangian of the GN model in Euclidean space is 

N 2 / ^ \ 2 

-Cgn = ipj + mo) V'i - y ( XI 

with 7 o = (Ti and 71 = (T 3 implying 7 /^ = 7 ^ = 7 ^ and { 7 ^, 7 i/} = 2i5^i^. Chiral symmetry is realized in a discrete way, 
V’i 75The term proportional to mo breaks chiral symmetry explicitly (it is analogous to a quark mass term). 
Therefore, in this work it is always set to zero, mo = 0 (similar choices are also implemented for the other models 
studied by us). 

Spontaneous symmetry breaking is only realized in the limit iV —> 00 [l^ . since for any finite N spontaneous 
symmetry breaking is excluded in 1 -b 1 dimensions . The chiral condensate arises upon condensation of the 

scalar-isoscalar field combination 'i/'jV't’ ^ ~ (where a sum over j is implied). 

In the limit N —>■ 00 analytic solutions for thermodynamical quantities including inhomogeneous condensation have 
been found @ (see also the discussion in Sec. uni). 


xGN model 

A straightforward extension of the GN model is obtained by adding a pseudoscalar term. The Lagrangian of the 
yGN model is 

N 2 

-^xgn = X - y 

This model contains a scalar field combination which corresponds to a tr-like particle, and a pseudoscalar field 

combination tpj^ which corresponds to an ? 7 -like particle. It is invariant under continuous CAi(l) chiral symmetry 
transformations, ipj —>■ and has certain similarities to one-flavor QCD (when the chiral anomaly is excluded). 


X 




( 2 ) 
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The xGN model is particularly interesting, because both the scalar and the pseudoscalar field configurations con¬ 
dense, when the temperature exceeds a critical value. The ground state is then a CDW 


NJL 2 model 

A further extension of the GN model is obtained by considering, in addition to the scalar-isoscalar field combination, 
three pion-like field combinations. In this respect the model is similar to two-flavor QCD. The Lagrangian of this 
so-called NJL 2 model is 


2 N 2 V / ^ ^ 

•^NJLa = EE ./ - y EE ) + ( E 

f=lj=l /=lLVj = i / \j=i 


(3) 


where / = 1, 2 is the flavor index. The model is invariant under chiral symmetry transformations SUl{2) x SUii{2), 

tpj,L -t ULlpj,L , -t Uuljjj^R , (4) 

with 

_ f V'i.i 


V'j = 




, tpj,L = Phi’j , i’j.R = Pr-iIjj , 


(5) 


and Pl and Pr are projectors onto left- and right-handed components, respectively. 

In contrast to the yGN model the ground state is not a CDW. Using the finite-mode approach we find that the 
phase diagram coincides with that of the GN model (cf. Sec. IVT). 


B. 3 + 1 dimensions: the NJL model 


The NJL model in 3 + I dimensions is one of the most famous effective chiral approaches to QCD. It has been 
extensively used in the vacuum and at non-zero temperature and density to study the spontaneous breaking of chiral 
symmetry and its restoration [cf. e.g. Refs. The Lagrangian (in the chiral limit for N colors and two flavors) 

is @,[13 


N 


^NJL = EE EKE ) + ( E 


/=ii=i 


N 


N 


/=i ■- ^i=i 


N 


1=1 


( 6 ) 


Chiral symmetry is realized in the same way as in the NJL 2 model [cf. Eqs. (|3]) and (O]. 

In the vacuum, the quark field obtains an effective mass, if the coupling constant G exceeds a critical value. 


6G 


N 


K E > 0 • 


(7) 


/=ii=i 


This effective mass is proportional to the chiral condensate in the vacuum, i.e., (Tq = —{N/6\/2G)m s ili3 , where 
the chiral condensate is defined according to 


N 


/=! 7 = 1 


( 8 ) 


In other words, the field combination which gives rise to a non-zero condensate is again When restricting 

this condensate to be constant, chiral symmetry restoration at high densities occurs via a first-order phase transition 
[ 71 I - I 73 . However, when allowing for an inhomogeneous condensate, the latter occurs at slightly smaller chemical 
potentials than the first-order phase transition. This is in agreement with the extended Linear Sigma Model results 
of Ref. H^j. 

Gontrary to the 1 + I dimensional models of Sec. Ill A[ the NJL model is not renormalizable. The equation for a 
takes the form 


tTo = —2\/2A^TOg I with I = i 


d'^p 


d?p 


1 


(27r)"‘ _ ^*2 


(27r)3 


(9) 
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where the integral I corresponds to a closed quark loop, i.e., to a tadpole diagram arising from the quartic NJL interac¬ 
tion of Eg. ( 1^ . which affects the quark propagator at the resummed one-loop level in the Hartree-Fock approximation 
[see Refs. [a-[l3 ^ detailed derivation]. 

The integral I is, however, quadratically divergent. Indeed, the NJL model is properly defined only after a regular¬ 
ization scheme has been chosen and a corresponding high-energy scale enters as a new parameter. Strictly speaking, 
each choice of regularization corresponds to a different version of the NJL model. Once the regularization has been 
fixed, the quantity mg in Eq. © and, as a consequence, all the relevant thermodynamical quantities, are finite. 

Especially in studies of the NJL model at nonzero density it is common to implement a three-dimensional cutoff 
[see e.g. Ref. [Ill] according to which 


I ^ Ik 


d\p\ f 

Jo 27r2 + TOg^ 


( 10 ) 


Note that the use of a four-dimensional covariant cutoff is possible for studies of the vacuum , but it is not easy 
to implement at nonzero temperatures and densities. 

However, a three-dimensional cutoff strongly suppresses the appearance of inhomogeneous phases. Namely, in an 
inhomogeneous phase such as a CDW the quark propagator is not diagonal and the ingoing and outgoing momenta 
can differ by a full wavelength. This hardly takes place when the momentum |p[ is limited by the cutoff A l7al. Hence, 
in order to realize a CDW, other regularization approaches must be used, such as the Pauli-Villars scheme or the 
proper-time regularization scheme SElzS- 

In this work we use the Pauli-Villars approach, which is a Lorentz (and gauge) invariant regularization procedure 
• It amounts to introduce Npy additional fictitious heavy fermions with mass Mk in such a way that the tadpole 
integral of Eq. m is modified according to 


/ —>■ /pv 


l^d\p\ I' 1 ^ ^ Ck ^ 


with the masses given by 


Idik — ^0^ + OfeApY , 


( 11 ) 


( 12 ) 


where Apy is the so-called Pauli-Villars high-energy scale. The constants Ck and ak are real dimensionless numbers, 
which are chosen in such a way that Ipv is finite. Let us show this explicitly for the case Vpy = 2. For large values 
of ff' the quantity in parentheses in Eq. (HU can be approximated by a Taylor expansion, 


(...) 


Then, by requiring 



1 + Cl + C 2 - 


1 ^0^(1^2) + Apy(aiCi -l- < 3 . 202 ) 

2 W 


o{\p\-^) 


(13) 


1 + Cl + C 2 — 0 and 0 . 1 C 1 H- 0 . 2 C 2 — 0 (14) 

the integrand of Eq. m falls off as \p\ ^ and is, therefore, convergent (although it explicitly depends on the scale 
Apy). Once Ipy is finite, the quark condensate, the quark mass, as well as all other relevant quantities are also finite. 
The conditions in Eq. da are met for ai = 2 and 02 = 1 with Ci = 1 and C 2 = — 2 . 

The procedure can be easily generalized to an arbitrary number of heavy fermions A^pvi 


A^pv iVpv 

1 + ^ Cfe = 0 , Y. • 


(15) 


For the case Vpy = 3 the previous equations are fulfilled by ai = 1, 02 = 2, 03 = 3 and Ci = —3, C 2 = 3 , 
C 3 = — 1. In Sec. IVII we will compute the phase diagram of the NJL model with inhomogeneous condensation using 
the Pauli-Villars regularization with two and three heavy fermions. 


III. FINITE-MODE REGULARIZATION OF THE 1 + 1 DIMENSIONAL GN MODEL 

In the following we discuss the finite-mode approach in detail, in particular its technical aspects, in the context of 
the 1-1-1 dimensional GN model in the large-V limit [cf. also Refs. m 0, liiill]. We reproduce the analytically 
known phase diagram, which exhibits an inhomogeneous crystal phase. 
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A. Partition function and Euclidean action 


The partition function of the 1+1 dimensional GN model m in Euclidean spacetime is 


Z = 


N 


n DiiJjD'iiJj 




,-Se [dj+j] 


with the action 


r N 


SEli’j.i’i] = / <Px 




i=i 


r 

2 


N 

1=1 


(16) 


(17) 


where ^ is the chemical potential. One can get rid of the four-fermion term by introducing a real scalar field cr, 


N 


z = / T'cr f n Dll’j Dll;j 


exp 


— / d‘^x 


2 ^^ 


N 

1=1 


(18) 


with the Dirac operator 


Q = l^id^ + 70/1 + cr . (19) 

Performing the integration over the fermionic fields results in 

J Da (detQ)^ exp ^ — . (20) 

Since detQ is real [s^l, (detQ)^ = det{Q^Q) > 0. Consequently, for even N 

Z = J , SE,efiW] = nJ ^' 2 : ^ ln[det(g+Q)]| , (21) 

where A = Ng^. Due to numerical reasons discussed in detail in Ref. [d^, when using the finite-mode approach it is 
highly advantageous to regularize the effective action expressed in terms of det(g^g) instead of the mathematically 
equivalent expression containing detQ. 


B. Finite-mode regularization, homogeneous condensate a — constant 


For numerical calculations it is convenient to work exclusively with dimensionless quantities. Therefore, we express 
all dimensionful quantities in units of ao which is the non-vanishing value of the constant condensate a at temperature 
T = 0 and chemical potential g — 0. The resulting dimensionless quantities are denoted by a hat e.g. 

T = T/ao, fi = /r/cTo, d = aja^, etc. 

We consider a finite spacetime volume with temporal extension Lq = L^ao (corresponding to the inverse temperature 
1/T = Lq) and spatial extension Li = Liao. The fermionic fields are expressed as superpositions of plane waves with 
periodic boundary conditions in spatial direction and antiperiodic boundary conditions in temporal direction. 


? /- . , lpj{xo,Xi) 

lfjj{Xo,Xi) = 

l/fXQ 

7 - N tpj{xo,Xi) 

lj’j{xo,Xi) = — 

v^o 

with discrete momenta 


^—i{kQXo+kiXi) 

(22) 

/ , Vj,no,ni p, —^ 5 

no,ni V 

^+r{kQXo+k-i_xi) 


^ . Vj,no,ni 1 -^—^ 5 

no,ni V 

(23) 


^ ZTT 

ki = —ni , no, ni G N , 

Li 



where r 7 j,no,ni and r]j,no,ni are dimensionless Grassmann variables. 
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For a homogeneous condensate a = constant, ln[det((5^Q)] can be expressed as a product over the modes {ko, fci) 
ln[det(( 5 'l'( 5 )] = In det - 70 /i + '^) ( “ + 7 oA + o’) 


ko,ki 


= In 


n [(*0 +kf + a‘^ - + (2/ifco)" 


ko,ki 


(24) 


Considering only a finite number of modes no = —Nq + 1, —Nq + 2,..., Nq — 1, Nq and ni = —Ni, —Ni + 1,..., A^i — 
1 , i.e., introducing momentum cutoffs 


27r 


2tt 


fccut ^ ^ 


1 


(25) 


Lq Li 

(chosen to be tt/Lq,! larger than the largest momenta considered) yields the finite-mode regularized effective action, 
SE,es{^) LoLia'^ 


N 


2A 

No Ni 


-) E E >" 

no = —N’o + l ni ——Ni 

27T^No{Ni + l/2)a'^ 

Nq Ni C r / 

- E E !■'{ N 

nQ — lni — — Ni ^L\ 


27T 

Lq 


no 


1 2 


27T 


-ni 




27T 

2 /r— no - 


1 2 


no - 1/2 

No 




Ni + 1/2 


I ~2 ~2 

+ cr - fi 


2/tfcj 


cut^o 1/2 


No 


(26) 


which is suitable for numerical evaluation. Minimizing this effective action with respect to a for various /t and T 
yields a{il,T), i.e., the “homogeneous phase diagram of the GN model” [6nl |. 

Infinite-volume continuum results are obtained in the limit Li —>■ 00 , —>• 00 , and fc™* —>■ 00 , which implies 

an infinite number of modes. Numerically one is, of course, restricted to a finite number of modes. In the following 
we discuss how to determine the parameters A, No, Ni, fc™*, and [Lq and Li are then given by Eq. ([25|l ] in 
an optimal way, i.e., how to obtain numerical results with a finite and rather limited number of modes, which are 
nevertheless very close to infinite-volume continuum results. 


C. Choosing and determining suitable parameters A, No, Ni, and 


The condensate a minimizes the effective action ((2fill . For ^ = 0 it is the solution of 

2tt^No{Ni + 1/2) 


0 = 4^^^= 2d 
da N 


\ Lcut Leu 


No Ni 

- 2 E E 

no — l ni =—ATi 


r,cut ''^o 1 /^ 
° No 


+ fci 


ni 


Ni + 1/2 


-b<T“ 


• (27) 


An obvious solution is d = 0. It corresponds to a minimum for T > Tc and to a maximum for T < Tc- For the latter 
case there are two additional solutions (corresponding to minima), which can be obtained from 


Q _ 27r^fVo(fVi -I- 1/2) ^ 2 yA 

\£.cut£.cut E—^ 

^^0 no = lni = -Ari 


Icut^o 1/2 

° No 




ni 


Ni + 1/2 


1 -1 


(28) 


To appropriately determine the parameters A, A^o, Ni, and fc™*, we consider and relate computations at /i = 0 
and a low temperature T Ri 0, where cr{T) r; ao, and at /x = 0 and the critical temperature T = To, where a just 
vanishes, i.e., a{T — e) > 0 and <j{T) = 0. The parameters A, A^i, and fc™* are the same for both simulations, 
while No = Noo for T m 0 and No = Nqc ^ Nqo for T = Tc- 

The parameters Noc, Noo, Ni, and can be chosen independently. The maximum number of modes oc NoqNi 
is, of course, limited by the available computer resources. Strategies for choosing these four parameters in an optimal 
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way, i.e., where systematic errors due to the finite spatial extension and the finite number of modes are minimized, 
are discussed in Secs. IIII C 11 IIII C'2l and IIII (JSl below. 

In contrast to that, and A cannot be chosen independently: = 2'kNqcTc [which follows from Eq. (1^51) ]. 

i.e., is related to iVgc. Since Tc is a priori unknown, setting fc™* to an appropriate value is a non-trivial task. 
Similarly, A depends on fVoc, iVi, and via Eq. ([28|) at T = Tc, where tt = 0, 


2/Ar|-|/r)A 1 Nqc Ni r / 2 / \ 2-| —1 

TT (iVi + 1/2) _ ^ ^ [ / z^cut^O - 1/2 


IVoc 

U 1 nQ — lni— — Ni ^ 


N, 


Oc 


kl'^ 


ni 


Ni + 1/2 


(29) 


To determine (without knowing Tc), we consider Eq. (l28l) also for T Ri 0, where ct ~ 1, i.e.. 


Nqo Ni 


+ 1 / 2 ) _ 1 

U i nQ — ln\— — N\ 


>ut ^o - 1/2 

fVoo 


'cut 


Ni + 1/2 


n -1 


+ 1 


(30) 


Since the left-hand sides of Eqs. (l29l) and (l30l) are identical, we can equate their right-hand sides and eliminate A, 


E 


ni— — Ni 




no - 1/2 

Nc 


2 



»1 \ 

Ni + 1/2/ 


2l 


-1 


1 

No 


E 


"0 = 1 



no - 1/2 

-^00 


2 



\ 2 

ni \ 

iVi + 1/2 j 


1 


= 0 . 


(31) 


Eor given Nqc, No, N, and one has to solve this equation to obtain fc™*. Then, A can be calculated using either 
Eq. (|2^ or (|3U1) . 


1. Optimizing Noc 


The numerically obtained critical temperature Tc = k™^/2 ttNoc should be insensitive with respect to variations 
of Nc, when keeping the other parameters fixed, in particular No- The corresponding optimal iVg/’* is, therefore, 
defined as the value of Nc which minimizes 



(32) 


(since Nqc S N, the derivative d/dNoc has to be understood as a finite difference). 

To study this optimization of Nqc independently of any error due to the finite spatial momentum cutoff and 
the finite spatial extension Li = 2 'k{Ni l/2)/kf^^, we consider for a moment the limit —>■ oo and Li —>■ oo 

(implying —>■ oo). In this limit 




ni 


—i k\ 




dki 


N 


Ni 

E 




N + l/2 

1 ni — — Ni 


N + 1/2 ’ fVi + 1/2 

Inserting these relations into Eq. (l311l and solving the integral results in 


iVi + 1/2 

£.cut 


dki . 


(33) 


1 

Nc 


Noc 

E 

"0 = 1 


tcut 

ftg 


no-1/2 A 

Nc ) 


-1 


1 

Noo 



no - 1/2 

No 


2 - 1 - 1/2 


+ 1 


= 0 . 


(34) 


As previously Eq. (IHTl) this equation has to be solved to obtain fc™*, which now only depends on Noc and No- 
In FIG. [T] we study the corresponding T/ = fc™*/27r/Vgc as a function of Nc (left panel) and Nc/No (right panel) 
for Noo e {64,128,256}: 

• For sufficiently large Noo and a suitably chosen Noc the resulting Tc should be close to the analytically known 
infinite-volume continuum result T/ = c^/tt r; 0.566 (where C denotes Euler’s constant) [^. One can clearly 
see that there are plateau-like regions, where this is the case. 


































• For a small number of temporal modes Nqc there are strong deviations, because the temporal momentum cutoff 
is rather small, = 2ttNocTc (for small Nqc, curves obtained with different Nqq fall on top of each other when 
plotted versus Nqc oc fc™*)- 

• For Nqc/Noo^O-^ there are also strong deviations, because the temperature corresponding to A^oo temporal 
modes, Tq = TcNqc/Nqq, is a poor approximation of zero temperature (for Nqc/Noq>0.2, curves obtained with 
different Nqq fall on top of each other when plotted versus Nqc/Noo = Tq/Tc). 

fc 

0.572 r 

0.571 - 
0 . 570 ;. 

0 . 569 ; 

0 . 568 ; 

0 567 ‘. . . . . 

0.566 ; 

0.565 . 

0 10 20 30 40 50 60 

FIG. 1: Tc as a function of Noc (left panel) and Noc/Nqo (right panel) for Nqo = 64 (green dots), A^oo = 128 (red dots), and 
Noo = 256 (blue dots) obtained from Eq. (I34II (the black lines indicate the infinite-volume continuum result Tc = e'^/rr). 

In other words, to obtain accurate results, 1 ^ Nqc <C iVoo has to be fulfilled, which is only possible, if a sufficiently 
large number of temporal modes A^oo are used. According to the definition (|32ll the optimal Nqc for given Nqq is the 
minimum of the corresponding curves in FIG. [1] 

In TABLElJone can see the accuracy I — Tc/ (c^/tt) of the numerically obtained Tc, for various values of Aqq. Note 
that even with a comparatively small number of Nqq = 32 temporal modes, the error is less than 0.1%. Also listed 
in TABLE U are A'g^*, the corresponding temporal momentum cutoff = 2ttNq^^Tc, and the corresponding 

temporal extension Lqq* = Nqq/TcNq^^ approximating T = 0. When increasing Nqq, there is a similar increase in 
the temporal momentum cutoff but only a slight increase in the temporal extension Aqo*- This is typical for 

lattice calculations, where cutoff effects are only polynomially suppressed, while finite-volume effects are exponentially 
suppressed. 


0 


T, 

0.572 
0.571 
0.570 
0.569 
0.568 
0.567 
0.566 
Noc 0.565 




0.05 0.10 


0.15 


0.20 0.25 


0.30 


035Noo 


Noo 

1 - fc/U 


Nof 

rcut,opt 

Kq 

0.0 
0 o 

32 

-6.320 • 

lO"'* 

6 

2.139 • 

10* 

9.401 

64 

-1.981 ■ 


11 

3.919 • 

10* 

10.26 

128 

-6.158 • 


20 

7.125 • 

10* 

11.29 

256 

-1.881 • 

lO"^’ 

36 

1.282 • 

10^ 

12.54 

512 

-5.650 • 

10”*’ 

66 

2.350 • 

10^ 

13.68 

1024 

-1.672 • 


120 

4.275 • 

ItF 

15.05 

2048 

-4.886 • 

10"' 

222 

7.908 • 

Tip' 

16.27 

4096 

-1.413 • 

10"' 

412 

1.468 • 

llF 

17.53 


TABLE I: Accuracy of the numerically obtained Tc for various values of Noo- 


2. Optimizing 


In 

(|M 1 ) 


the following we investigate the error associated with a finite spatial momentum cutoff fc™*, i.e., instead of Eq. 
we return to Eq. and solve this equation to obtain fc™* for given Nqc, Nqq, Ni, and fc™*. Similarly to 
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Sec. IIII CD we define the optimal 


as the value of which minimizes 



(35) 


for given Nq^, Nqq, and A^i. 

We choose Nqq = 256, the corresponding optimal Nqc = 36 (cf. TABLE |l| and various numbers of spatial modes 
A^i. In FIG. Owe study T), = l:™*/27r7Voc as a function of (left panel) and fc™*/(7Vi + 1/2) (right panel) obtained 
from Eq. dSI]) for £ {64,128, 256}: 


fe fc 



FIG. 2: Tc as a function of (left panel) and /{Ni + 1/2) (right panel) for Ni = 64 (green curve), Ni = 128 (red curve), 
and Ni = 256 (blue curve) obtained from (1311) with Nqq = 256 and Nqc = 36 (the black lines indicate the infinite-volume 
continuum result Tf, = e'^/vr). 


• For sufficiently large Ni and a suitably chosen fc™* the resulting should be close to the analytically known 
infinite-volume continuum result Tc = c^/tt Ri 0.566. Again one can observe plateau-like regions, where this is 
the case. 

• For a small spatial momentum cutoff fc™* there are strong deviations (for small fc™*, curves obtained with 
different A^i fall on top of each other, when plotted versus fc™*). 

• For fc™*/(A^i -I- 1/2) > 1.0 there are also strong deviations, because the extent of the periodic spatial dimension 
Li = 2'k{Ni + 1 /2) /fc™* is quite small and, therefore, a poor approximation for infinitely extended space 
(for fc™*/(Afi -b 1/2) >1.0, curves obtained with different A^i fall on top of each other, when plotted versus 
fcr*/(iVi + 1/2) oc 1/Li). 

In other words, to obtain accurate results, 1 <C fc™* A^i has to be fulfilled, which is only possible, if a sufficiently 
large number of spatial modes Ni is used. According to the definition (1551) the optimal fc™* for given Ni is the 
maximum of the corresponding curves in FIG. [2] 

In TABLE |IT] one can see the accuracy 1 — ///(c^/tt) of the numerically obtained Tc for various values of A^i, again 
for A^oo = 256 and the corresponding optimal A^oc = 36. Note that errors due to the finite-mode regularization in 
temporal direction and in spatial direction have opposite sign (cf. FIG. [T] and FIG. [5]). Consequently, one obtains the 
most accurate result for Tc not for Ni —> oo, but when both errors almost cancel each other. This is the case for 
Ni = 249, i.e., for Ni similar to Nqq (for a more detailed discussion cf. Sec. IIII C 31) . Also listed in TABLE n are 
the spatial momentum cutoff fc™*’°P* and the corresponding spatial extension = 2tt{Ni -\- l/2)/fc™*’°*’*. Again, 
when increasing A^i, there is a similar increase in the spatial momentum cutoff fc}'^*’°^*, but only a slight increase in 
the spatial extension approximating infinite volume. As already mentioned this is typical for lattice calculations, 
where cutoff effects are only polynomially suppressed, while finite-volume effects are exponentially suppressed. 
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Ni 

1 - 

^CUtjOpt 

i-opt 

32 

-h5.789 • 10-* 

2.273 - 10^ 

8.99 

64 

+1.766- 10-'‘ 

3.991 - 10" 

10.15 

128 

+4.324 • lO"'" 

7.138 - 10" 

11.31 

249 

-2.720 ■ 10"* 

1.261 - 10^ 

12.43 

256 

-1.003- lO"” 

1.291 - 10^ 

12.48 

512 

-1.460 - lO"*’ 

2.329 - 10^ 

13.83 

1024 

-1.798- lO"'^ 

4.179 - 10^ 

15.41 

2048 

-1.871 - lO"'" 

7.616 - 10^ 

16.90 

4096 

-1.882 - lO"*’ 

1.403 - 10^ 

18.35 


TABLE II: Accuracy of the numerically obtained Tc for Aqo = 256, A^oo = 36, and various values of A^i. 


3. Optimizing the ratio of Nqo and Ni 


As already mentioned the maximum number of modes oc A^oo-^i is limited by the available computer resources. In 
the previous subsection it has been observed that for Nqq = 256 and = 249 the errors in Tc due to the finite-mode 
regularization almost cancel. Also for other choices of Nqq such a nearly perfect cancellation is present for r; Nqq, 
as collected in TABLE Hill Moreover, note that the temporal momentum cutoff and the spatial momentum 

cutoff are close to each other. 

Similarly, in FIG. [3] we compare the numerically obtained Tc for various Nqq and Ni with the infinite-volume 
continuum result Tc = eP j-n. The figure suggests choosing Nqq = as a simple rule, which obviously leads to very 
accurate numerical results (the black filled circles in FIG. [3]). Unless mentioned otherwise, we will use Nqq = Ni in 
the following. Of course, such a cancellation of errors might not occur for quantities other than Tc- Moreover, when 
considering also the possibility of an inhomogeneous condensate, as we will do in Sec. IIIID1 it could be necessary to 
have a finer resolution or a larger extent of the spatial dimension, which might require a rather large A^i ^ Nqq. 


A'oo 

Ni 

^opt 

r cut,opt 

Kq 

^CUtjOpt 

32 

30 

6 

2.137- 

10^ 

2.190 - 

10^ 

64 

62 

11 

3.919 - 

10^ 

3.932 - 

10^ 

128 

125 

20 

7.124 - 

10^ 

7.034 - 

10^ 

256 

249 

36 

1.282 - 

10^ 

1.261 - 

10^ 

512 

497 

66 

2.351 - 

10^ 

2.287 - 

10^ 

1024 

994 

121 

4.310 - 

'W 

4.188 - 

IfF 

2048 

1985 

223 

7.844 - 

ItF 

7.710 - 

1(F 

4096 

3966 

414 

1.475 - 

liF 

1.428 - 

IfF 


TABLE III: Pairs Nqq and Ni, where the errors due to the finite-mode regularization for Tf, almost cancel, and the corresponding 
N°f, and 


4- Summary 

Based on these investigations we propose and adopt the following strategy to determine the parameters Nqc, Ni, 
and A: 

(1) Use Nqq = Ni as large as possible (limited by the available computer resources). 

(2) Determine the corresponding optimal Nqc from a computation in the limit A^i —>■ 00 and /c™* —> 00 as described 
in Sec. IIII CTl fcf. also TABLE HI third column). This computation also provides a value for (TABLE HI 
fifth column), which is a good approximation for fc™* at finite, but large Ni and fc™*. Assign this value to 

(cf. Sec. IIII G 31) . 
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T 

^ c 



FIG. 3: Tc as a function of A^i for several Noo (filled black circles indicate symmetric choices A^oo = A'^i, the black line is the 
infinite-volume continuum result Tc = 6^1%). 


(3) Solve Eq. (1511) to determine fc™* (now for finite Ni and A™*) for the previously chosen Nqq and Ni [step (1)] 
and Nqc and fc™* [step (2)]. 

(4) Determine A via Eq. (|2^ [or equivalently via Eq. (l30l) ]. 

For all further computations, e.g. when computing the phase diagram for a homogeneous condensate a = constant 
or an inhomogeneous condensate a = cr(a;i), the parameters fc™*, fc™*, and A are not changed anymore. The 
temperature T = 1/Lq = fc™‘/27rA^o can be adjusted by using different numbers of temporal modes Nq. Nq^ 
corresponds to the critical temperature Tc- 


D. Computation of the phase diagram for a homogeneous condensate a = constant 


To determine the phase diagram for a homogeneous condensate, i.e., cr as a function of the chemical potential /r 
and the temperature T, one proceeds as in the previous subsection. From 


0 = 


da N 

one can derive the generalization of Eq. (E51) for arbitrary fi> 0, 


(36) 


0 = 


27r^Afo(A^i-H/2) 

\ £.cut l.cut 

/AfbQ 


No Ni 

- 2 E E 

no = l ni— — Ni 




£.cut rap —1/2 1 . I ;„cut tii \ i 

0 No / I 1 Ni+1/2 j 


— p? -b I 2/lfc™* 


^ 0 - 1/2 
0 No 


■ (37) 


If this equation has a solution d'^ = A > 0 ior given {p,T), and if SE,es{'/A) < SE,es{0), this solution is the value 
of the chiral condensate, i.e., d — dz\/A, and iP,T) is inside the chirally broken phase. If there is no such solution, 
or if S'£;_eff('\/3) > SE,esiQ), then d = 0 and {p,T) is a point inside the chirally symmetric phase or on the phase 
boundary. 

The phase diagram for a homogeneous condensate obtained with Nqq = Ni = 192 (with corresponding Nqc = 28, 
fc™* = 9.974 • 10^, = 1.011 • 10^ and A = 0.3328) is shown in FIG. ID It is in excellent agreement with the 

infinite-volume continuum result 
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f 



FIG. 4: Phase diagram of the GN model for a homogeneous condensate cr = constant for Noo = Ni = 192 [the light gray cnrve 
is the infinite-volnme continuum result [^1. 


E. Finite-mode regularization, spatially inhomogeneous condensate cr = 

To study the possibility of a spatially inhomogeneous condensate, a is written as a superposition of a finite number 
of plane waves, 


&{xi) 


crjxi) 

CTo 


M 


E 

m— — M 


^-tpxi 



~ _ Stt - f \* 

P ^ 771 , C-\-rri — [C—jn) 

Li 


(38) 


as done in Eqs. (HH) and (1251) for the fermionic fields t/jj ■ The resolution of a should be coarser than the resolution of 
i’j, i.e., M Ni, to obtain stable and meaningful numerical results [cf. Ref. [i^ for a detailed discussion]. 

In the case of a spatially inhomogeneous condensate plane waves are no longer eigenfunctions of the Dirac operator 
Q. Consequently, ln[det((5’^(5)] cannot be expressed as a product over modes as done in Eq. (|25)) . One has to represent 
Qp as a matrix, where the rows and columns correspond to the plane-wave basis functions of the fermionic fields 
^^i(koxo+kixi)^^ [cf. Eqs. (1221) and (1251) ]. 


{ko,k\Q^Q\k'o,k[) = 


2 pl 

^ ^ / dxo I 

LqLi Jo Jo 


+M 


dxi e+<'^°^o+k^x,) _ ^ ^ ^ 


m— — M 




+M 


X ( +7m5/.+7oA+ E 

m' — — M 


-tp'xi \ 

Cm' — I C 

vE / 




+M 


. -f i7ofco + *7ifci+7 oA+ E 

•'0 m=-M V Ll 


+M 


X ( - *70^0 - *71^1 + 7oA + 


^-ip Xi \ 


E ) 

V-^1 '' 


, — ik'^xi 


(39) 


with p = 21 : 111 /Ll and p' = 27rm'/Li. These matrix elements can be calculated analytically. Note that this matrix 
representation of det((5’f(5) has a block-diagonal structure with 2No blocks of size 2(2fVi -|- 1) (the blocks are labeled 
by fco = ^oi tbe rows and columns of each block correspond to the spatial momenta ki and and the two spin 
components). ln[det(( 5 ^( 5 )] is then the sum over ln[det(...)] of the blocks, where each term of that sum can be 
evaluated numerically e.g. by means of an LU decomposition. 
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F. Computation of the phase diagram for an inhomogeneous condensate a — a{xi) 

When allowing for a spatially inhomogeneous condensate a = (t(xi), there are three phases: 

(I) For small chemical potential p, and low temperature T chiral symmetry is broken by a homogeneous condensate 
a = constant ^ 0 [corresponding to cq ^ 0, Cm = 0 for m 0 in Eq. (l38l) ]. 

(II) For high temperature T chiral symmetry is intact and a = 0 [corresponding to Cm = 0 for all m in Eq. (1381) ]. 

(Ill) For large chemical potential p and low temperature T there is a spatially inhomogeneous condensate a = a{xi) 
[corresponding to 7 ^ 0 for at least one m 7 ^ 0 in Eq. (1551) ]. 

The phase diagram for a spatially inhomogeneous condensate obtained with Nqq = = 192 and M = 10 (with 

corresponding iVoc = 28, fc™* = 9.974 • 10^, = 1.011 • 10^ and A = 0.3328) is shown in FIG. [S] It is in excellent 

agreement with the infinite-volume continuum result . 

f 



FIG. 5: Phase diagram of the GN model for a spatially inho mog eneous condensate a = CT(a;i) for Noo = W = 192 and M = 10 
[the light gray curve is the infinite-volume continuum result [48.1^1. 

The numerical determination of the phase boundaries is discussed in detail in Ref. (d^ and, therefore, only sum¬ 
marized briefly in the following. 

• Phase boundary I-II: 

The phase boundary between a = constant 7 ^ 0 (phase I) and <7 = 0 (phase II) can be determined as explained 
in Sec. IIII Dl for the phase diagram for a homogeneous condensate. 

• Phase boundary I-III: 

To determine the phase boundary between cr = constant ^ 0 (phase I) and the inhomogeneous crystal phase 
(phase III), one again has to find the minimum of S'e^eff/A/ with respect to <7 as a function of (/i, T). This time, 
however, <7 is not a constant, but a superposition of plane waves [cf. Eq. (1381) ]. The minimization has to be done 
with respect to the coefficients c„. 

• Phase boundary II-III: 

To determine the phase boundary between (7 = 0 (phase II) and the inhomogeneous crystal phase (phase III), 
one can in principle proceed as for the phase boundary I-III. Note, however, that inside the crystal phase in 
the vicinity of the phase boundary I-III the constant <7 of the phase diagram for a homogeneous condensate is a 
local minimum (the corresponding phase transition is of first order), while in the vicinity of the phase boundary 
II-III it is a saddle point (the corresponding phase transition is of second order). Therefore, a computationally 
simpler and cheaper way to determine the phase boundary II-III is to study the smallest eigenvalue of the 
Hessian matrix 

) C2m — 110 ( 077 . 1 ) 7 C'2m-\-l — Icn(Cy 77 ) , (^ 0 ) 

(7—0 

with 0 < m,m' < M. A negative eigenvalue amounts to a direction of negative curvature and, therefore, 
indicates the existence of an inhomogeneous condensate. 




d d SE,es{^) 


dCrr,. dC„ 


N 











14 


Since the finite-mode approach allows to determine the condensate a at given temperature T for arbitrary chemical 
potential it is straightforward to study and reproduce the order of the transition along the phase boundaries I-II, 
I-III, and II-III. 


IV. PHASE DIAGRAM OF THE 1 + 1 DIMENSIONAL yGN MODEL 

We proceed in the same way as explained in detail in the previous section for the GN model. After introducing two 
real scalar fields a and t], the partition function of the xGN model can be written as 

Z = J Da y a; cr^-f 77 ^^ - i ln[det(Q'''Q)]| , (41) 


where 


Q = 'Ifj.dfj. -I- 7oAi + cr + mi5- 


We then apply the finite-mode regularization, i.e., in analogy to Eq. (1551) the scalar fields a and rj (the condensates) 
are represented as a sum over a finite number of modes, 


a{xi) 


Vixi) 


a{xi) 

So 

So 


M 

■n— — M 
M 

h=-M 






x/Ta 


e+m — (C—m) 


d-|-m — (d—771) , p — - TTl 

Li 


(42) 

(43) 


Similary, Q^Q is written as a matrix, where the rows and columns correspond to plane-wave basis functions 
^^i(koxo+kixi) LqLi, [cf. Sec, nil El and Eqs. (l22l) and (l23l) for details]. Since Se. eff is invariant under the transfor¬ 
mation (cr, 7 ) —>■ R{a,r]) with R G 0(2), dimensionful quantities are expressed in units of Eg, where 



T=0,M=0 



(44) 


and denoted by a hat '. 

We have studied the phase diagram of the xGN model using M = 10 modes for the condensates and Nqq = Ni = 96 
modes for the fermionic determinant. 

For temperatures T > Tc = /tt and arbitrary chemical potential /i chiral symmetry is restored, i.e., the effective 

action SB) is minimized for Cj = dj = 0, which corresponds to vanishing condensates a = rj = 0. 

For T < Tc we find several local minima of SE,eS, which are given by Cm = iidm ^ 0 for a single mode m, while 
Cj = dj = 0 for all other modes, i.e., j ^ m. The condensates a and rj are harmonic functions, i.e., CDWs, with the 
same amplitude, but with a relative phase shift ± 7 r/ 2 , implying E = constant [cf. Eq. (1441) ]. The minimal values of 
SE,eS are plotted in FIG. [51 as functions of the chemical potential p for m = 0,1, 2,3 and two different temperatures 
T = 0.378,0.189 (Vg = 24,48, while Vgg = Vi = 96 and M = 10). For ^ = 0 the absolute minimum of «S'£;,eff 
corresponds to m = 0, i.e., a,rj = const. The wavelength is proportional to /r, i.e., for increasing p the absolute 
minimum of SE,eS corresponds to larger and larger m > 0. Two examples of the resulting CDWs (T = 0.095 < Tc 
[Vg = 96] and fi = 0.295, 0.875) are shown in FIG. jT] 

The resulting phase diagram is, therefore, quite different from the phase diagram of the GN model: there are only 
two phases, for T < Tc a, CDW, and for T > Tc chiral symmetry is restored. This is in agreement with analytically 
known results 


V. PHASE DIAGRAM OF THE NJL 2 MODEL 


Again the technical steps needed to compute the phase diagram closely parallel those discussed for the GN model 
and the xGN model. This time four real scalar fields a and tt^, j = 1,2,3 are required, where 

^ = J ( n , S'E,eff[cr, TTj] =N J I -I- ^ i ln[det(Q'l'Q)] | , 


(45) 
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SE.eft SE,et{ 



FIG. 6: xGN model. Local minimal values of Se.bS as functions of the chemical potential jj, for m = 0 (blue), m = 1 (red), 
m = 2 (green) and m = 3 (gray) and two different temperatures T — 0.378 {Nq = 24, left panel) and T = 0.189 {No = 48, right 
panel). 


&,fj &,fj 



FIG. 7: xGN model. Resulting chiral CDWs for T = 0.095 < Tc {No = 96) and jX = 0.295 (left panel) and ji = 0.875 (right 
panel). 


with 


3 

Q = +7oAi + o-+ *75 ■ (46) 

1=1 

cr and TTj as well as det((5^<3) are then finite-mode regularized as done in Eqs. (l38l) . (l42l) . and (l43l) and Sec. IIIIEl re¬ 
spectively. Due to the invariance of SE,eS with respect to (cr, tti , 7r2, tts ) —>■ R{a, tti , 7r2, tts) with i? G O(4), dimensionful 
quantities are expressed in units of 



T= 0 ,/i =0 



(47) 


and denoted by a hat ~. 

We have studied the phase diagram of the NJL 2 model using M = 4 modes for the condensates and iVoo = 
Ni = 72 modes for the fermionic determinant. For any temperature and chemical potential the four condensates are 
proportional to each other, i.e., cr oc tti oc 7r2 oc 713 oc E and also proportional to the chiral condensate cr of the GN 
model. Consequently, we obtain exactly the same phase diagram for the NJL 2 model as for the GN model, which is 
shown in FIG. [5] These findings extend existing results, where only a CDW [cr(a;) oc cos(26a:) and 713(2;) oc sin(26a;)] 
has been considered [s^. Our results show that an inhomogeneous phase is indeed present at larger ^ and not too 
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large T. However, this inhomogeneous phase exhibits solitonic structures and not a CDW. This result is similar to 
those found for the NJL model in 3 + 1 dimensions, as discussed e.g. in Sec. rvTl or Ref. . 


VI. PHASE DIAGRAM OF THE NJL MODEL 

We investigate the phase diagram of the NJL model at nonzero temperature and density under the assumption that 
only the chiral condensate a defined in Eq. ([ 8 ]) condenses. Thus, we do not take into account condensation of the 
pion-like field combinations appearing e.g. in Eq. ([5]). The chiral condensate a is, in general, a function 

of the three spatial coordinates x = (xi,X 2 ,x^), i.e., a = a(x). However, previous investigations based on special 
ansatze have shown that modulations in more than one dimension are not favoured energetically [3^ [^. Thus, for 
the sake of simplicity in the following we assume that a depends only on one of the three spatial coordinates, e.g. X 3 , 
i.e.. 


(T = a{xz) 


o-Q 


M 


E 

m— — M 



C+m — (C—m) 5 P — 


27T 

^3 


(48) 


Of course, a study of the NJL model in the context of the finite-mode approach without this assumption is an 
interesting topic which we plan to investigate in the future. 

Proceeding as in Sec. lHI i^ for the GN model we obtain the partition function of the NJL model in 3-1-1 dimensions, 

Z = J , SeM<^] = nJ “ ln[det(QtQ)]| , (49) 

where Q = Jfj.dfj, J- 70 /^ + xtiq. The Pauli-Villars regularization can be implemented by adding heavy fermions as 
explained in Sec. IHB[ 

P ATpv 

<5'£;,eff,Pv[o'] = S'_E,eff[o'] -N d'^x ^ Ck ln[det((5j)Qfc)] , (50) 

k=l 


with Qk = 7^9^ J- 7 oAi + Mk- 

For a convenient comparison with the existing literature on the NJL model (and in contrast to the previous 
sections, where we discussed U-l dimensional models) we express our results in the following in units of GeV. To 
this end, the two parameters of the NJL model, the coupling constant G and the Pauli-Villars energy scale Apy, are 
fixed by requiring that a certain effective mass tuq is realized (we perform computations for three different choices, 
mg G {250 MeV, 300 MeV, 350 MeV}) and that the pion decay constant reproduces the correct value in the chiral 
limit, = 88 MeV [Sll, [H^. [For the evaluation of the pion decay constant in the framework of the NJL model, we 
refer to Ref. 0.] 

The resulting phase diagrams for Npy = 2 additional heavy fermions, effective quark masses TOq G 

{250MeV, 300MeV, 350MeV} and M = 5 and Nqq = Ni = 120 modes are shown in FIG. The corresponding 
Pauli-Villars cutoffs are Apy = {736.8 MeV, 647.4 MeV, 608.7 MeVj, the maximum momentum used in the expansion 
of the chiral condensate is approximately (M/V)fc{'^* Ri 650 MeV, where fc™* denotes the maximum momentum in 
the expansion of the fermionic fields, see Eq. (I25|) . 

Let us discuss the role of the two energy scales Apy and /c™* in somewhat greater detail. As mentioned already in 
Sec. H.B, the NJL model is non-renormalizable. As a consequence, the NJL model is only defined once the regular¬ 
ization has been fixed. Moreover, the corresponding energy scale (the Pauli-Villars cutoff Apy in our implementation 
of the NJL model) should be regarded as a physical parameter of the NJL model, with crucial impact on all physical 
quantities. On the other hand, our numerical approach contains also a cutoff due to the lattice regularization in 
momentum space, i.e., due to using a finite number of modes, see Eq. (j25p . The cutoff fc™* is purely technical and 
should be taken as large as possible. In particular, should be much larger than Apy. In the numerical calculations 
leading to Fig. 8 we have that fc™*/Apy = 25.4, thus is about 16 GeV. This is indeed a very large value that 
assures that for all practical purposes the results do not depend on Such a high value also assures that the 

continent of Fig. 8 is not an artifact of our numerical calculation. We have also verified that relevant quantities, such 
as the critical temperature, do not depend on once it is chosen sufficiently large. 

Quite remarkably, these phase diagrams are similar to that obtained for the 1-f 1 dimensional GN model (cf. FIG. 

This result suggests that the 3-1-1 dimensional NJL model, which represents a non-renormalizable, but in many 
aspects realistic chiral model of QCD, generates a phase diagram whose most salient features can be understood in a 
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simpler 1+1 dimensional field theory. However, it should also be stressed that this result is obtained in the specific 
case of the Pauli-Villars regularization and that the existence of an inhomogeneous condensate depends on the chosen 
regularization scheme. 


T [GeV] T [GeV] 



FIG. 8: Phase diagram of the NIL model for a spatially inhomogeneous condensate a = (t{xz) for A^pv = 2, M = 5, 
Nqq = Ni = 120 and three different effective quark masses mg = 250 (red), mg = 300 (green) and mg = 350MeV (black). 
There are three phases: (I) a homogeneons chirally broken phase o = constant ^ 0, (II) a chirally restored phase a — 0, (III) 
an inhomogeneons phase (the left panel is a zoomed version of the right panel). 

Independent of the concrete choice of the effective quark mass ttiq there is an inhomogeneous phase for large chemical 
potential fi and small temperature T, termed “continent” in Ref. [s^. However, at smaller /i the detailed shape of 
the phase diagram depends on the value of the effective quark mass. While at small TOq an inhomogeneous “island” 
may be separated from the “continent”, at larger mg the island and the continent merge. Note that our results in the 
region /i < 0.4 GeV are in agreement with the findings of Ref. [3l|. At somewhat larger /i they agree with the results 
of Ref. [1^, although the outlines of the continent were not traced to very large /r in that work. We find that, at even 
larger values of /r, the transition temperature between the chirally restored and the inhomogeneous phase decreases 
with /i, which is similar to the GN model. 

We have repeated the study of the phase diagram of the NJL model for different volumes and verified that the results 
are very stable upon changing the volume. In particular, the form of the continent remains practically unchanged 
when Nqq is modified. We have checked the cases Nqq = Ni = 48, 72, 96,120 and already for Nqq = Ni = 48 the 
curves look just as in Figs. 8 and 9. We have then selected for the plots the highest number of modes {Nqq = Ni = 120) 
with which we could perform our numerical study in a reasonable computational time. 

We have also compared regularizations using Npy = 2 and A^pv = 3 additional heavy fermions (cf. FIG. IH). The 
effect of Npy on the shape of the phase diagram is rather mild: at intermediate chemical potential /i and temperature 
T the inhomogeneous continent becomes somewhat larger when using a larger number of regulators. 

Although a detailed study of the order of phase transitions at the phase boundaries is in principle straightforward, it 
would require a substantial numerical effort, especially in the case of a second-order phase transition, as expected for 
instance at the H-HI boundary. We have, however, verified that smoothly approaches zero at the H-HI boundary. 
On the other hand, for small and constant temperature, is a slowly decreasing function of the chemical potential 

/i. 


VII. SUMMARY AND OUTLOOK 

In this work we have used the finite-mode approach to investigate numerically the emergence of inhomogeneous 
chiral condensation in effective quark models in 1 +1 and 3 + 1 dimensions. The main aim has been the determination 
of inhomogeneous condensation in QGD-inspired models without using a specific ansatz. 

We have shown that our method accurately reproduces well-known analytical results concerning the phase diagram 
and the inhomogeneous condensation of 1 + 1 dimensional models, in particular the Gross-Neveu and chiral Gross- 
Neveu model. By applying the approach to the NJL model in 3 + 1 dimensions we could reproduce previous results 
based on specific ansatze for small chemical potential. In addition to that we were able to show that the inhomogeneous 
continent of Ref. extends to very high densities, but not to arbitrarily large temperatures. Due to the fact that 
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FIG. 9: Phase diagram of the NJL model for a spatially inhomogeneous condensate for A^pv = 2 (blue) and Npv = 3 (black) 
(further parameters chosen as in FIG.[H]|. 


our results for the NJL phase diagram differ from previous ones at high density, it would be an interesting task for 
future studies to confirm or to falsify the presence of the continent using different approaches. 

It is interesting to note that the GN model and the NJL 2 model have the same phase diagram, which in turn is 
also very similar to the NJL model in 1+3 dimensions. On the contrary, the phase diagram of the yGN model is 
completely different, showing that a different flavor structure has an impact on the phase diagram. 

Our approach following Refs, is based on expanding fermionic fields and condensates in terms of plane waves. 

In this respect it is quite similar to existing general methods to compute inhomogeneous condensates [s^, . 

In other aspects it is, however, quite different. For example our method is based on a minimization of the effective 
action, not the Hamiltonian. A detailed and direct comparison regarding the computational efficiency of our method 
and those discussed in Refs. [13, [13,113 is difficult and would require computations of the same quantities within the 
same model using both types of methods. Since the basic idea of expanding in terms of plane waves is the same, 
we expect both types of methods to perform on a quantitatively similar level. Therefore, further studies would be 
necessary to clarify why the shape of the “continent” (see Fig. 0) looks different in the two methods. 

Note that within our approach it is straightforward to replace the plane waves by another set of basis functions, 
as discussed and numerically demonstrated in Refs. [H, Ull. A first and straightforward step would be to still use a 
plane-wave expansion for the condensate, but to focus exclusively on higher modes in regions where the condensate 
is expected to exhibit strong oscillations (typically in inhomogeneous regions at large chemical potential /i). Using a 
small number of modes with suitably chosen wave number might allow to explore the phase diagram quite extensively 
at rather moderate computational cost. The choice of the wave numbers could even be automated, i.e., included in 
the minimization procedure for the effective action. Another possibility within our numerical framework is to study 
specific ansatze to identify a small set of highly relevant degrees of freedom. This could not only provide certain 
physical insights, but also reduce the computing time significantly. The latter might be of particular importance, 
when studying two- or three-dimensional variations of the condensates. 

In the future one could also apply the finite-mode approach to study the phase diagram of purely hadronic the¬ 
ories such as the extended Linear Sigma Model Si- This approach is capable to correctly describe the vacuum 
phenomenology as well as the nuclear matter ground state properties [sJl. As already shown in Ref. [l^l, an inho¬ 
mogeneous condensate in the form of a chiral density wave is favored with respect to a constant condensate at high 
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density (> 2po, where po is the nuclear saturation density). It is an open question whether other structures minimize 
the effective potential even further. More generally, one could also apply the finite-mode approach to quark-based 
sigma models, e.g. Refs. [H, [s^. 

Other interesting projects are the study of higher-dimensional modulations beyond the ansatz used in Ref. [s^ 
and, since there is also no limitation in the number of inhomogeneous fields in the finite-mode approach, the study 
of interweaving chiral spirals [s^. Further effects at high densities such as inhomogeneous diquark condensation in 
1 -I-1 as well as 3 -I-1 dimensions can also be taken into account . Moreover, the models that we have studied were 
investigated in the chiral limit only. Future work could thus include a non-zero bare quark mass into the effective 
approaches. 
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